Method and apparatus for modelling a tire for use with a vehicle spindle-coupled simulator

ABSTRACT

A method for developing a tire model for use with an effective road profile in testing an automotive vehicle on a spindle-coupled road simulator defines a flat surface road plane in a three coordinate reference system to represent the effective road profile with the three coordinates comprising a plane vertical deflection, a plane radial contact angle, and a plane lateral contact angle. The tire is then excited over a predetermined range while on a tire test apparatus for developing the tire model. The tire test apparatus includes a flat tire contact plane moveable so as to contact the tire, the flat tire contact plane defining a three coordinate reference system to represent the effective road profile in the three coordinates.

FIELD OF THE INVENTION

The present invention relates to tire models for automotive vehicle simulators, and particularly for spindle-coupled vehicle simulators.

BACKGROUND OF THE INVENTION

The use of large-deflection tire models to predict vehicle loads for design or simulation testing has become an important application in the automotive vehicle development process. A wide range of tire models have been used to simulate the tire. In general, current tire models can be divided into three categories: simple models, ring elastic foundation (REF) models, and detailed models.

Two simple tire models are shown in FIGS. 1 and 2. These models assume that the tire envelopes an obstacle and that the effect the tire tread has on contact and deformation is negligible. In a radial spring model (FIG. 1), the tire envelopes the terrain with springs that are radially attached to center of the tire. Another simple model, commonly referred to as an ADAMS model (FIG. 3), uses an "equivalent ground plane" to calculate the longitudinal and vertical spindle forces. In both models, the net spindle force is generated by summing the individual spring forces in the vertical and horizontal directions.

In contrast to simple models, a Ring Elastic Foundation (REF) model represents the tire as an elastic foundation supporting a tread. REF models can be represented either with partial differential equations or with finite elements. Detailed tire models are three dimensional finite element representations of the complete tire (FIG. 4). These models are more representative of an actual tire than the previously discussed models, but a major disadvantage of this type of model is the large number of degrees-of-freedom needed and consequently the intensive computing time required.

Many existing tire models use a finite element analysis (FEA) approach to determine vehicle loads. However, these currently available models have inaccurately estimated tire loads through computer simulation. In addition, all of the tire finite-element models discussed above are time-intensive to build and run. Vehicle spindle-coupled simulation requires a tire model that is relatively easy to model from readily-measurable parameters, which is not computationally intensive, and which is easy to implement with a spindle-coupled simulation controller utilizing an effective road profile.

SUMMARY OF THE INVENTION

The present invention overcomes the disadvantages of the related art by providing a method and apparatus for modeling a vehicle tire using system identification techniques for use with an effective road profile. To develop such a tire model for an effective road profile driven vehicle simulation process using system identification techniques, a new method and apparatus for acquiring input and output data, and new conventions for defining an effective road profile are required. Since models developed using system identification methods are particularly useful for relating output data to input data (or vice versa if the model is inverted), a coordinate system is selected to represent the road surface that make it generic from a durability testing standpoint while also allowing it to be physically introduced using an experimental test system. A `flat surface road plane` coordinate system is defined to satisfy both of these requirements. The flat surface road plane coordinate system is defined in three coordinates: a plane vertical deflection at the tire patch center (ERP₋₋ 1), a plane radial contact angle with the tire (ERP₋₋ 2), and a plane lateral contact angle with the tire (ERP₋₋ 3). Loads which affect the structural durability of vehicles can thus be reproduced at the test tire by achieving a correct spindle force and acceleration in the vertical and longitudinal axes of a tire, in addition to a correct spindle force-moment relationship in the lateral axis of a tire. Experimental data for use with the tire model system identification technique is therefore obtained in the form of spindle force and acceleration, and spindle force-moment, collectively the outputs, as well as the input tire plane vertical deflection, plane radial contact angle, and plane lateral contact angle.

A tire test stand capable of producing the required input and collecting the desired output with reference to the flat surface road plane coordinate system is also necessary. The tire test stand generates random flat surface radial deflections and flat surface lateral contact angles with a test tire mounted thereon. After the radial deflection input/spindle force output data is generated, it is used to develop a radial tire model using the system identification technique.

Thus, an advantage of the present invention is a method for modeling a vehicle tire for use with a spindle-coupled simulator which is relatively easy to model from readily-measurable parameters, which is not computationally intensive, and which is easy to implement with a spindle-coupled simulation controller utilizing an effective road profile.

Another advantage is a tire test stand which is capable of providing test inputs to a vehicle tire and generating test outputs in a flat surface road plane coordinate system required for to develop a tire model for use with an effective road profile.

BRIEF DESCRIPTION OF THE DRAWINGS

These and other objects, advantages, and features of the present invention will be apparent to those skilled in the vehicle testing arts upon reading the following description with reference to the accompanying drawings, in which:

FIG. 1 shows a prior art simple tire model;

FIG. 2 shows a second prior art simple tire model;

FIG. 3 shows a third prior art simple tire model, referred to as an ADAMS model;

FIG. 4 shows a finite element representation of a tire for use with a finite element tire model;

FIG. 5 shows an effective road profile for use in combination with a tire model to eliminate sample road testing of a vehicle to be tested on a spindle-coupled simulator;

FIG. 6 is a schematic flow diagram of the effective road profile control method for use with a spindle-coupled road simulator;

FIG. 7 is a schematic flow diagram of an effective road profile control method for a spindle-coupled simulator according the present invention;

FIG. 8 shows development of an effective road profile generic to a vehicle from a set of unique vehicle dependent responses and a tire model;

FIG. 9A and 9B show, respectively, a top view and a side view of a tire test stand according to the present invention for development of a tire model for use in spindle-coupled road simulator;

FIG. 10A and 10B, respectively, show a front view and side view of a tire on a flat tire contact plane of a tire test stand in accordance with the present invention which varies in radial displacement, radial contact angle, and lateral contact angle;

FIG. 11 and 12 show, respectively, a tire freebody force and parameter model and a tire model diagram;

FIG. 13 is a first schematic flow diagram showing development of a tire model through use of a system identification technique for use with an effective road profile control method;

FIG. 14A and 14B is a second schematic flow diagram showing a second effective road profile control method for a spindle-coupled road simulator; and

FIG. 15 is a third schematic flow diagram showing a third effective road profile control method for a spindle-coupled simulator.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

Referring now to the drawings and in particular to FIG. 5 thereof, a vehicle A has a suspension configuration which includes a spring 70 above a spindle 72. A spindle force history 73a and spindle acceleration history 73b are generated when vehicle A is driven over a test road, such as is typically available at a vehicle proving grounds. The spindle force and acceleration histories 73a, 73b are used in conjunction with a model 74 of a tire attached to vehicle A when driven over the test road to develop an effective road profile 75. For purposes of this disclosure, an effective road profile is a set of signals which, in conjunction with a tire model, and when used as a feedback parameter in a standard response simulation test process, provides accurate excitation of the spindle-coupled simulator to simulate a test road. A vehicle B, having a different configuration, for example an air suspension 76, can be accurately tested without the need for driving the vehicle B over a test road. That test is accomplished by using the effective road profile 75 and a second tire model 77, which represents a tire used on the vehicle B, to generate a spindle force and acceleration histories 78a and 78b which would be developed had vehicle B been driven over the test road with a tire model 77. Thus, only one vehicle configuration need be driven over a test road to obtain the data necessary to develop an effective road profile which can subsequently be used for testing other vehicle configurations.

A schematic flow diagram of the just described process is seen in FIG. 6. In ellipse 82, vehicle dependent parameters, such as spindle force and spindle acceleration, are obtained for a first, or test, vehicle when driven over a test road. The vehicle dependent parameters are then used in conjunction with a tire model 83 to yield a vehicle independent effective road profile 84. Since the effective road profile is independent, that is, it is representative of the actual test road without regard to vehicle configuration or tire model, it may be used with a model of a spindle-coupled simulator 85 in a feedback control process to develop a simulator drive file 86 which can be used as an input to a spindle-coupled simulator 87 with a second vehicle coupled thereto. It is important to note that the second vehicle need not have been driven over the test road.

Turning now to FIG. 7, a more detailed flow diagram of a preferred embodiment of the effective road profile control method for a spindle-coupled simulator of the present invention is shown. Beginning at the top left of FIG. 7, a proving grounds facility 91 has a test road over which a vehicle A, indicated by box 92, is driven. Spindle dynamics data from the vehicle A is collected and is represented at cylinder 93, the data including spindle force and spindle acceleration as previously discussed. Typically, a vehicle has four spindles and spindle dynamics data is collected from each spindle. The data may be collected using force and acceleration transducers, or from other data collection devices known to those skilled in the art.

The collected dynamics data is then taken to a test facility or laboratory, as indicated by the dotted line 94, where it is used in conjunction with an inverse of a tire model, indicated at 95, to estimate a generic effective road profile, indicated at 96. The inverse tire model 95 is based upon a tire model of a tire used on vehicle A when driven over the test road. Those skilled in the art will recognize that the steps previously described are accomplished for each spindle and tire of the vehicle A, and that a single tire model for vehicle A may be used for each spindle.

The tire model, discussed in more detail below, can be represented as a transfer function developed from a set of linear differential equations which represent a physical tire. The model must be invertible, that is, have an inverse as that term is known by those skilled in the systems dynamics and control arts. The inverse of the tire model is needed for use with the spindle dynamics data since it is the tire input which must be predicted, although the tire model was developed to produce the spindle response to a known input.

Still referring to FIG. 7, the effective road profile 96 is used in a feedback control loop as indicated by line 97. The feedback control loop includes a summer 98 at which signals from the effective road profile 96 are compared to an output, referred to as "remote" response and indicated generally at 100. The remote response 100 is an output from a spindle-coupled road simulator transfer function HBT, indicated generally at 102, which includes an inverse of a tire model 104 of a tire on a second, or sample, vehicle B, indicated at 106. Development of the tire model for a tire of vehicle B is accomplished similarly to that for the tire of vehicle A, as described below. It should be understood that development of the tire models, both for vehicle A and vehicle B, need not be done at any particular time or in any particular order but only that such models be available when required by other steps of the method of the present invention.

Vehicle B, which has a different configuration than that of vehicle A as described above with reference to FIG. 5, is coupled to a spindle-coupled road simulator 108. A servo-controller 110, such an a conventional servo-control apparatus, is used for control of the road simulator (FIG. 7). System modeling of the vehicle B, the road simulator 108, and the servo-control 110, collectively comprise a transfer function HB, indicated generally at 112. Spindle dynamics data for vehicle B, indicated at cylinder 114, are generated from the road simulator 108 with vehicle B attached thereto. The spindle dynamics data 114 is then used in conjunction with the inverse tire model 104 for vehicle B to yield the remote response 100. The transfer function HB, the spindle dynamics 114, and the inverse tire model 104 collectively comprise the transfer function HBT.

As previously discussed, a sum of the generic effective road profile 96 and the remote response 100 is calculated at 116 and is operated on by an inverse transfer function HBT⁻¹, indicated generally at 118, before being written to a drive file B indicated generally at cylinder 120. The set of control signals thus generated in the drive file B represents a unique control input for a particular simulator-vehicle-controller combination. The vehicle B can then be further tested while coupled to the spindle-coupled simulator without the need for driving it over the test road.

Turning now to FIG. 8, it is seen that a tire model transforms vehicle dependent spindle responses, which are unique to a particular vehicle, into an effective road profile, which is vehicle independent. Various tire models could be used to simulate the tire providing that the particular model used has an inverse. Simple tire models typically assume that the tire envelopes an obstacle and the effect the tire tread has on contact and deformation is negligible. A net spindle force is generated by summing individual spring forces in the vertical and horizontal direction. For example, a radial spring model envelopes the terrain with springs that are radially attached to the center of a tire.

Ring elastic foundation models represent the tire as an the elastic foundation supporting a tread. These models are represented either with partial differential equations or with finite elements. This model requires fewer degrees-of-freedom than other finite element models to produce accurate results.

Detailed tire models are three dimensional finite element representations of a complete tire. These models are more representative of an actual tire and therefore, predict loads more accurately.

Regardless of the level of detail, all of the finite-element model forms discussed above are time-intensive to build and run and require an understanding of the tire material properties. If such requirements are met, the previously discussed models must be invertible, or so constructed to provide the requisite bridge between vehicle dependent spindle responses and an effective road profile. Preferably, however, a tire model which is easily derivable from readily measurable parameters, and which is easily implemented with the spindle-coupled simulator controller is desirable.

To develop such a tire model, system identification methods are preferably used to identify dynamic properties of a vehicle tire to be modeled. The system identification modeling approach requires experimental input/output test data which allows the dynamics of a system to be observed, and aids in the selection of an appropriate model form. Alternatively, a model form can be assumed, based upon an understanding of the known physical characteristics of the system. Once a model form for the system is selected, the model parameters (e.g., coefficients of the selected differential equations) can be identified based upon a system identification method (such as least squares estimation) that minimizes the output error between the physical system and modeled system for a given input.

This input/output test data is obtained through use of a tire test stand such as that seen in FIGS. 9A and 9B. The tire test stand has a three point, fixed reacted load frame with a shaft connected thereto supported by zero-lash tapered roller bearings to prevent tire wind-up. A wheelpan is used to provide excitation to the tire and serves as a flat contact plane which varies in radial displacement, radial contact angle, and lateral contact angle (FIGS. 10A and 10B). The system identification technique thus uses the tire contact plane, or wheelpan, as a coordinate system in which to develop a set of dynamically varying data including radial deflection, radial contact angle, and lateral contact angle of a tire to be modeled.

Referring again to FIGS. 9A and 9B, a tire test stand for use in developing a tire model for a vehicle to be coupled to the spindle-coupled simulator is shown. A tire 130 is mounted on a spindle 132 which is journalled in spindle support bearings 134 on opposite sides of the tire 130. A set of supports 136 attach the spindle to a tire support structure 138. A pair of double axis, support force transducers 140 are mounted on the tire supports 136 for collecting spindle dynamics data. The tire 130 is in contact with a tire contact plane 142, modeled as a flat surface road plane. Lateral contact angle input is provided through a piston 144 connected to the tire contact plane 142, the piston 144 preferably being hydraulically actuated. Radial deflection is provided by a piston 146, also connected to the tire contact plane 142. The piston 146 is also preferably hydraulically actuated.

The tire test stand includes measurement transducers for measuring wheelpan displacement, wheelpan force and acceleration, tire patch to wheelpan displacement, and spindle force. Additionally, a displacement transducer, preferably a laser displacement transducer, can be mounted in the wheelpan to track the tire patch trajectory in order to develop a loss of contact profile.

The tire test stand has a hydraulic actuator with a capacity in the range of 15 KIPS, a 3 inch displacement, and a three stage servo-valve with a flow rate of approximately 90 gallons per minute. The actuator/valve peak velocity is in the range of 100 to 150 inches per second. Resonance frequencies for the tire test stand are preferably higher than the spindle-coupled simulator controller bandwith. Example frequencies for the tire test stand are 70 to 100 hertz for the wheelpan lateral mode, and 80 to 100 hertz for the wheelpan vertical mode. Those skilled in the art will recognize that tire test stand system capacities, displacements, and valve flow rates and peak velocities can vary according to application and that the values given are by way of example and not by way of limitation.

Turning again to FIGS. 10A and 10B, a coordinate system is selected to represent the road surface so as to make it generic from a testing standpoint while also allowing it to be physically introduced using an experimental test system. As such, a "flat surface road plane" coordinate system is defined to satisfy both of those requirements. The flat surface road plane coordinate system is defined in three coordinates: a plane vertical deflection at the tire path center (ERP₋₋ 1), a plane radial contact angle with the tire (ERP₋₋ 2), and a plane lateral contact angle with the tire (ERP₋₋ 3). Loads which affect the structural durability of vehicles can thus be reproduced at the test tire by achieving a correct spindle force and acceleration in the vertical and longitudinal axes of a tire, in addition to a correct spindle force-moment relationship in the lateral axis of a tire. Experimental data for use with the tire model system identification technique is therefore obtained in the form of spindle force and acceleration, and spindle force-moment, collectively referred to as the outputs, as well as the input tire plane vertical deflection, plane radial contact angle, and plane lateral contact angle.

A system model can be derived using physical principles and observed experimental behavior, as shown in FIGS. 11 and 12. In FIG. 11, a free body force diagram is shown for a tire. Discussion of the tire model is best understood with reference to FIG. 12, and expanded tire diagram. A spindle force, FSP, acts on the spindle-coupled tire net mass, M, which is derived from physical principles, at the tire spindle, indicated at 132 in FIG. 11. A non-physical modal tire model matches experimental observation and is represented as modal mass MTM with a first spring 140 mounted between the net mass M and the modal mass MTM, and also having a spring 141 between the modal mass MTM and the wheelpan 144. Both springs 140, 141 have a spring stiffness constant K. Likewise, a damper 142 is modeled in parallel with the spring 140 between the tire modal mass MTM and the net mass M, and a damper 143 is in parallel with the spring 141 between the tire modal mass MTM and the wheelpan 144. Both dampers 142, 143 have a damping constant of C. Displacement parameters of the spindle, YSP, the modal tire, YTM, and an effective road profile, YERP, are shown in FIG. 12.

A Laplace transfer function representing the tire system model is:

    M.sub.tm s.sup.2 F(s)+2Cs F(s)+2K F(s)=C.sup.2 s.sup.2 Y(s)+2CKs Y(s)+K.sup.2 Y(s)

where s is a complex variable used in the Laplace transform to reduce differential equation representations of dynamical systems to simple algebraic form. This transfer function can be converted from the continuous Laplace s-domain form to the discrete difference equation z-domain form for use with computational methods, as on a computer, by substitution using the backward-difference relationship:

    s=(1-z.sup.-1)/T

where T is the discrete computational time interval.

S=(l-z⁻¹)/T

K=spring stiffness constant,

C=damping constant,

F(s)=Spindle force,

Y(s)=Tire deflection, and

M_(tm) =tire system mass

This transfer function can be converted to difference equations by using the backward difference relationship for s, as defined above, to convert from a continuous to a discrete form for use with computational methods, as on a computer. As those skilled in the art will recognize, many techniques are available to then estimate the parameters of the resulting linear system model. A preferable method is a least squares algorithm.

Referring now to FIG. 13, a flow diagram for a system identification technique used to identify the dynamic properties of a test tire is shown. A broadband excitation 151 of the tire test stand 152 is required to accurately determine tire dynamic properties. In order to excite the tire over its full dynamic range without damage, the proper input excitation range must first be established. To establish this input excitation range, a tire is first mounted in the test stand and a static tire deflection necessary to achieve 1/4 of the vehicle weight is determined. Next, a maximum deflection amplitude is determined by moving the wheelpan into the tire until a predetermined maximum spindle force is achieved at a specific tire pressure. These ranges are used to determine maximum displacement which will be used to achieve desired force levels based on static load and deflection conditions.

Once the maximum deflection amplitude is determined, the tire test stand system is excited by a servo hydraulic actuator, which is computer controlled, to identify steady state behavior of the tire. Preferably, a frequency between 0 and 300 hertz of random white noise is used having an amplitude with a weighting function of the inverse of a square of the frequency. To avoid exceeding static tire deflection, the mean-centered white noise is preferably passed through a quadratic transformation. This allows the input displacement to reach desired levels while simultaneously insuring that tire patch compressive forces are maintained within the static and dynamic limits.

Since tire temperature has a significant effect on its apparent pressure, it is preferable to excite the tire with a small amplitude input, for example approximately a 1 inch signal, for a given time, for example, approximately 20 minutes, prior to exciting the tire on the test stand and taking the requisite measurements. Such a step allow tire pressure to stabilize thereby providing more accurate measurements.

Referring again to FIG. 13, broadband tire input signals 153 and spindle force response signals 154 are collected during the broadband excitation 151 of the tire test stand 152. Both broadband tire input X1 and spindle force response Y1 are then input to the system identification 155. The system identification 155 is a least squares parameter estimation to minimize the model output error against measured output for the model selected. Parameter 156 determinations, such as the spring stiffness constant K and the damping constant C, are subsequently output from the system identification 155 to the tire model 157. The tire model 157 also receives broadband tire input X2 and calculates a spindle force response 158 based upon the tire model. The model spindle force response 158 is then compared to the actual spindle force response Y2 from the test stand to provide tire model correlation 159.

If spindle response data, such as spindle acceleration and spindle force, are not directly obtainable when driving a vehicle over a test road, other vehicle response data can be recorded and used to derive spindle responses through an intermediate process, as shown in FIG. 14A. This intermediate process may be necessary when time and expense constraints prevent direct collection of the spindle dynamics data. For example, instrumentation of a test vehicle to obtain the spindle dynamics data may require time consuming vehicle disassembly and assembly to properly locate the appropriate transducers on or near the vehicle spindles. Beginning at the top left of FIG. 14A, a proving grounds facility 161 has a test road over which a vehicle A, indicated at box 162, is driven. The vehicle response data which, for example can be force and acceleration data from a component on the vehicle in the vicinity of the spindle, such as the chassis, is collected, as represented at cylinder 163. Preferably, such data is collected in the vicinity of each of a vehicle's four tires. As previously discussed, such data may be collected using force and acceleration transducers or other data collection devices known in the art.

The vehicle response data is then taken to a test facility or laboratory, as indicated by the dotted line 164, to derive spindle dynamics data, and ultimately an effective road profile. As a first step in that direction, the vehicle response data 163 is used in a feedback control loop, generally indicated at 165 to produce a drivefile A, indicated at cylinder 166, which is used to drive a spindle-coupled road simulator with Vehicle A, which is the vehicle driven over the test road at the proving grounds, coupled thereto. A servo-controller 167, such as a conventional servo-control apparatus is used for control of the road simulator (FIG. 14A). System modeling of the vehicle A, the road simulator 168, and the servo-control apparatus 167, collectively comprises a transfer function HA, indicated generally at 169. A "simulated" response, indicated generally at 170, is produced as a result of the road simulator and is compared at summer 171 to the vehicle response 163. The resulting difference between the vehicle response 163 and the simulated response 170 is used as a control signal after being operated on by an inverse transfer function HA-⁻¹, indicated generally at 172.

Continuing with FIG. 14B, the data developed for the drivefile A, at cylinder 166, is then used an input the spindle-coupled simulator and the transfer function 169 in order to collect spindle dynamics data, indicated at cylinder 173, which includes spindle force and spindle acceleration. This spindle dynamics data is then used in conjunction with an inverse of a higher model for vehicle A, indicated at 174, to estimate a generic effective road profile indicated at cylinder 175. As described above, the inverse tire model 174 is based on a tire model of a tire used on vehicle A when driven over the test road. Once a generic effective road profile has been developed, it can be used as described above with reference to FIG. 7. Briefly, a difference between the effective road profile 175 and a "remote" vehicle response 180 is calculated and operated on by an inverse transverse function HBT⁻¹, indicated generally at 176, before being written to a drivefile B, indicated generally at cylinder 177. Vehicle B, which has a different configuration than that of vehicle A, is coupled to the spindle-coupled road simulator 178. A servo-controller 179 is used for control of the road simulator 178. System modeling of the vehicle B, the road simulator 178, and the servo-controller, indicated generally at 179, collectively comprise a transfer function HB, indicated at 181. Spindle dynamics data for vehicle B, indicated at cylinder 182, are generated from the road simulator 178 with vehicle B attached thereto. The spindle dynamics data 182 is then used in conjunction with the inverse tire model, indicated at 183, to yield the remote vehicle response 180. Vehicle B can thus be tested, and road forces simulated, without having to drive the vehicle B over a test road.

Turning now to FIG. 15, a second alternative embodiment for the present invention is shown. As with the previously described embodiments, a vehicle A, indicated at 200, is driven over a test road at proving ground facilities, indicated at 201. Spindle dynamics data at 202 is collected and used with an inverse tire model, indicated at 203, to develop a generic effective road profile 204. Unlike the previously discussed embodiments, however, the effective road profile 204 is input directly to a controller of a spindle-coupled simulator, indicated generally at 220. The controller 220 has an inverse tire model 222, for a tire of vehicle B, embedded directly into the controller 220. Such a configuration allows the effective road profile to be used directly as a control signal for the controller 220, thus eliminating certain computational steps required by the previously described embodiments when vehicle B is coupled to the spindle-coupled simulator (see FIGS. 7 and 14A, 14B).

Although the preferred embodiments of the present invention have been disclosed, various changes and modifications may be made thereto by one skilled in the art without departing from the scope and spirit of the invention as set forth in the appended claims. 

We claim:
 1. A method for modeling a tire for use with an effective road profile in testing an automotive vehicle on a spindle-coupled road simulator, the method comprising the steps of:(a) defining a flat surface road plane in a three coordinate reference system having three coordinates to represent said effective road profile, said three coordinates comprising:a plane vertical deflection; a plane radial contact angle; and a plane lateral contact angle; (b) mounting a test tire in a tire test stand against a contact plane movable in each of said three coordinates, with said tire, said contact plane, and said test stand comprising a tire test system; (c) exciting said contact plane with said tire in contact therewith over a predetermined excitation range with respect to each of said three coordinates; (d) collecting a set of excitation input data and a set of response output data from said tire test system, said data referenced to said three coordinates; (d) matching said input data with said output data; (e) estimating a spring stiffness constant of said tire system from said input data and said output data based on a linear model of said tire system; (f) estimating a tire system mass and a tire system damping constant from said input data and said output data using a system identification technique based upon a second order model of said tire system and a spring constant of approximately twice said spring stiffness constant estimated in step (e) so that a tire model is developed, (g) exciting said tire test system with a predetermined input signal to generate an actual spindle force from said tire; (h) estimating a simulated spindle force from said tire model based upon said predetermined input signal; and (i) comparing said actual spindle force with said simulated spindle force to verify said tire model with respect to said effective road profile.
 2. A method according to claim 1 wherein said excitation range is determined by mounting said tire in said test stand and measuring a static tire deflection necessary to achieve 1/4 of the weight of a vehicle, and by moving said tire contact plane into said tire until a predetermined maximum spindle force is achieved at a predetermined tire pressure to determine a maximum deflection amplitude.
 3. A method according to claim 1 wherein a transfer function for modeling said tire system in step (f) is

    M.sub.tm s.sup.2 F(s)+2Cs F(s)+2K F(s)=C.sup.2 s.sup.2 Y(s)+2CKs Y(s)+K.sup.2 Y(s)

where, s=(l-z⁻¹)/T, K=spring stiffness constant, C=dampint constant, F(s)=Spindle force, Y(s)=Tire deflection, M_(tm) =Tire system mass, and z=a discreet difference domain variable.
 4. A method according to claim 1 wherein said tire test system is excited over said excitation range between a frequency approximately 0 and 300 hertz of random mean-centered white noise having an amplitude with a weighting function of the inverse of a square of the frequency.
 5. A method according to claim 4 wherein said mean-centered white noise is passed through a quadratic transformation to avoid exceeding static tire deflection.
 6. A tire test apparatus for developing a tire model for use with an effective road profile in testing an automotive vehicle on a spindle-coupled road simulator, the apparatus comprising:a support structure; a spindle attached to said support structure for receiving a tire thereon; a flat tire contact plane moveable so as to contact said tire, said flat tire contact plane defining a three coordinate reference system to represent said effective road profile, said three coordinates comprising:a plane vertical deflection; a plane radial contact angle; and a plane lateral contact angle; contact plane lateral adjustment means attached to said tire contact plane for varying a lateral contact angle of said tire contact plane with respect to said tire; contact plane radial adjustment means attached to said tire contact plane for varying a radial contact angle of said tire contact plane with respect to said tire; and transducer means for measuring contact plane displacement, contact plane force and acceleration, tire to contact plane displacement, and spindle force.
 7. A tire test apparatus according to claim 6 wherein the contact plane lateral adjustment means and the contact plane radial adjustment means comprise, respectively, a lateral adjustment hydraulically operated piston and a radial adjustment hydraulically operated piston.
 8. A tire test apparatus according to claim 6 wherein the transducer means for measuring said tire to contact plane displacement is a laser displacement transducer.
 9. A tire test apparatus according to claim 6 wherein the tire contact plane is a wheelpan. 